LIMITS ON NEUTRINO RADIATIVE DECAY FROM SN1987A 
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ABSTRACT 

We calculate limits on the properties of neutrinos using data from gamma- 
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■ ray detectors on the Pioneer Venus Orbiter and Solar Max Mission satellites. A 

. massive neutrino decaying in flight from the supernova would produce gamma 

rays detectable by these instruments. The lack of such a signal allows us to 

Q^l constrain the mass, radiative lifetime, and branching ratio to photons of a 

^ ' massive neutrino species produced in the supernova. 

1. Introduction 
• 1—1 , 

^ I The occurence of Supernova 1987a in the Large Megellanic Cloud has proven to 

' be among the most fruitful experiments in the heavenly laboratory for the confirmation 

of "known" physics and the constraining of new physics. Aside from its obvious impact 
upon the study of the late stages of stellar evolution in general and upon supernova 
physics in particular, models for SN87a have become an industry for the study of the 
couplings of light particles (neutrinos, axions) to ordinary matter. In this work, we 
discuss limits upon the properties of neutrinos independent of specific model for the 
supernova. 

When a supernova occurs, the bulk of the binding energy of the progenitor star 
(~ 3 X 10^^ erg) is released in neutrinos, a fact predicted by theory and confirmed by the 
observation of a neutrino burst from the supernova, with a characteristic temperature 
of about Tu ^ 4.5 MeV. If at least one species of neutrinos is unstable and if it couples 
to the photon, then some of these neutrinos will decay to photons en route, which 
are potentially detectable as MeV gamma rays. At the time of the supernova burst's 



arrival at earth and environs, there were several satellites operating in the solar system 
capable of detecting the decay photons in the course of their watch for gamma ray 
bursts. Analyses of the data from one of these detectQrs, on board the Solar Max 
Mission (SMM) satellite, has already been presented.cm Here, we examine the data 
from Gamma Burst Detector on the Pioneer Venus Orbiter (PVO). 

2. Expected Gamma- Ray Signal 

We assume that 1/3 of the total supernova energy is released in a species of 
massive neutrinos. These neutrinos then decay in flight into a photon and a light 
neutrino {e.g., an electron neutrino). Typically, the decay products will each have 
energies of 1/2 x 3 x T^, ^ 7 MeV. 

We will consider decays of the form u — * u'+'j. We expect that the parent neutrino 
will be a massive exotic neutrino, while the daughter neutrino will be a member of a 
light family such as Ug. In particular, this decay allows two possibilities: the helicity of 
the daughter neutrino may either be flipped or the same with respect to the parent as 
the photon takes away a unit of spin (assuming both neutrinos are relativistic) . From 
quantum mechanics, then, we know that the distribution of the photon in the rest 
frame of the parent will be proportional to either cos^(0/2) or sin^(0/2), where is 
the rest-frame angle between the directions of the parent neutrino and the photon. 

At a given time at the detector, photons are received that resulted from neutrinos 
decaying on the surface of an ellipsoid defined such that the sum of the time between 
the supernova and the neutrino decay and the time between the decay and the detection 
of the photon is equal to the time of detection. Actually, this is only approximately 
true — there is an additional delay incurred due to the finite mass of the parent neutrino. 
Furthermore, due to the relativistic beaming of the decay products into the forward 
direction, many more photons are received from long, skinny ellipsoids (corresonding 
to short delay times) than from others (corresponding to longer delays). Thus, the 
relativistic delay will play a correspondingly larger role in these events. 

From relativistic kinematics, the angular distribution of daughter photons in turn 
determines the distribution of photon energies as a function of neutrino energy, for each 
of the two possible decays (flip or no flip). The time delay is determined by the energy 
of the neutrino and the decay ellipsoid above. Thus, the photon spectrum as a function 
of time is determined. It is given by 

dN = *^ ' fi(E,fi)dtdtddEdfi 

AtxD^ 7r 



(1) 



where B^ is the branching ratio of the parent neutrino to photons, is the neutrino 
number flux at the detector, is the differential number luminosity as a function of 
E, the parent neutrino energy, D = 1.7 x lO^^cm is the distance to the supernova, td 
is the time of the neutrino's decay, v is the neutrino's velocity and 7 = E/my is the 
relativistic factor, r is the neutrino lifetime, t is the time of the photon's arrival at the 
detector, is the cosine of the "lab-frame" decay angle between the parent neutrino's 



direction and the photon, and fi{E,^) is the distribution of angles as a function of 
neutrino energy for each of the two hehcity possibihties. The Dirac 6 function enforces 
the arrival ellipsoid. In this paper, we will concentrate on the limit of short decay times, 
td ^ IT <^ D, so the delta function becomes 6 {t — td{l — vji)). 

Changing variables from /i to the photon energy /c, and integrating over the decay 
time td gives 

dN = 5^-^e-2'=*/'"''"^/i(E, k) dkdEdt, (2) 

where fi{E,k) now gives the distribution of photon energies as a function of parent 
neutrino energy: 

f(j^ f,) _hk/E^ no flip , . 

^^^'^'-\2{E-k)/E' flip. 

Similar results have been derived for related cases before.!^ Note that this distribution 
is a function only of the combination of neutrino parameters rriiyT and B^; to break this 
degeneracy between mj, and r, we must relax our assumption of short decay times. 

We will also assume that the initial neutrino luminosity is given by a zero- 
chemical-potential Fermi-Dirac distribution. Normalized to the known temperature and 
total energy of the suprenova, this gives 

120 E' 

LdE) - ^^Y^^- (4) 

Finally integrating over this distribution, this gives the photon spectrum as a 
function of time, 



dN _ 240 Et 2k 



e-^f'i/^^^hiik/n). (5) 



dkdt rn^AnD'^T^m^T 

where hi is a separate function for each of the helicity possibilities that is of order unity 
for the energies and temperatures of interest. 

3. PVO Instrumentation 

The Pioneer Venus Orbiter Gamma Burst Detector (OGBD)i was designed to 
detect gamma ray bursts — transient, high energy events that last from milliseconds to 
tens of seconds. It has four separate bands, 0.1-0.2 MeV, 0.2-0.5 MeV, 0.5-1 MeV, and 
1-2 MeV. In the background mode that we will be using, it has a timing resolution of 
12 or 16 seconds with full spectral information {i.e., counts in each of the four bands). 

Because gamma-ray bursts are singular events, the OGBD was designed to have 
full-sky coverage, although it is incapable of independently providing directional infor- 
mation about the photons it receives. However, the OBGD detectors do, in fact, face 
the South Ecliptic Pole, and thus are ideally suited for measurements of gamma-rays 
from the LMC. Thus, even though it was not "pointing" at the LMC at the time, it 
can still be used to analyze data associated with the supernova. Although the response 
of the instrument changes for is a function of the gamma-ray direction, the bulk of the 
decay photons come from the direction of the supernova itself on the sky (the excep- 
tions are the rare photons which reach us after decaying at a very large angle from the 
outgoing neutrino). 



4. Analysis 

After a brief look at the data from the four OGBD channels in Fig. 1, it is clear 
that there is no obvious signal of gamma rays over the background (the variation in 
the signal is in fact consisitent with pure noise). Thus, we will be putting limits 
on our free parameters: rriuT and B^. 



Fig. 1. Raw counts/sec in each of the four OGBD channels for the time surrounding the 
arrival at Venus of the supernova's optical burst, at time UT — 27325 = sec. 

To do this, we utilize the SPANAL code which was developed for analyzing the 
instantaneous spectra of gamma-ray bursts. Given the spectrum averaged over 12- or 
16-second bins, SPANAL calculates the best fit set of parameters, m^r and B^, for 
the data (thus far, we simulataneously fit the data for up to six time bins, or about 
90 seconds after the supernova). In Fig. 2, we show the results for a small area of 
parameter space. The points are minima, with ~ 22 (with 2 de grees of freedom) 
for all of the points. For the point m^r = 10^ keV sec, we show in addition the 1-a 
upper limit on B^; it is not very far from the minimum, so we expect that the 
actual upper limits are not far from these points throughout. The area below the line 
is roughly excluded; the area above is allowed by the OBGD data. To compare with 
previous analyses, Kolb and TurnerB found a limit in a similar area of parameter space 
of B^ < 2.8 X 10~^° without taking into account the angular distribution of decay 
photons (by simply assuming k = E/2). 



In order to more fully map the limits of parameter space, we must keep the full 
information of the arrival ellipsoid in the delta function of Eq. 1. For such long arrival 
times that this effect matters, however, the expected number of gamma rays is much 
lower (since fewer neutrinos would have decayed before reaching earth), so the limits 
will be somewhat weaker. We are currently pursuing the analysis in this realm. 



Fig. 2. minima for the fit of the theoretical spectrum of photons (Eq. 5) to the OGBD 
data after the supernova. For m,yT = 10^ keV sec, we also show the 1-a upper limit. 
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